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ABSTRACT 

It  is  well  known  that  the  boundary  layer  flow  over  a  surface  with  a  region  of  concave 
curvature  is  susceptible  to  centrifugal  instabilities  in  the  form  of  Gortler  vortices.  In  the  limit 
of  large  Gortler  number  (a  parameter  which  is  a  measure  of  the  curvature  of  the  surface) 
the  effect  of  a  crossflow  component  in  the  underlying  basic  flow  has  been  shown  to  stabilises 
these  modes  and  thus  render  the  Gortler  vortex  mechanism  inoperable  in  these  situations. 
Here  we  consider  the  effect  of  crossflow  when  the  Gortler  number  (and  the  scaled  spanwise 
wavenumber  of  the  vortex)  are  both  order  one  quantities.  The  parabolic  partial  differential 
equations  governing  the  linear  evolution  of  a  Gortler  vortex  in  a  three-dimensional  boundary 
layer  are  solved  numerically.  Our  results  suggest  that,  at  least  for  small  magnitude  crossflows, 
the  Gortler  vortex  instability  mechanism  is  still  operable.  In  addition  we  consider  the  effect 
of  an  applied  pressure  gradient  within  the  boundary  layer  on  the  instability  mechanism  and 
demonstrate  that  a  favourable  pressure  gradient  renders  the  boundary  layer  more  susceptible 
to  the  Gortler  vortex  instability;  this  is  in  stark  contrast  to  the  case  of  Tollmien-Schlichting 
waves  where  a  favourable  pressure  gradient  stabilises  the  flow. 
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§1  Introduction 

There  are  many  mechanisms  which  can  promote  transition  in  a  boundary  layer  and 
as  such  it  is  important  to  be  able  to  predict  how  a  given  physical  characteristic  will 
effect  the  breakdown  of  the  flow.  Our  concern  is  with  the  problem  of  how  two  such 
features  will  compete  to  change  the  transition  characteristics  of  the  flow.  The  two 
features  considered  here  Eire  crossflow  and  the  concave  curvature  of  the  plate. 

Many  physical  situations  have  flow  components  in  more  than  one  direction,  one 
example  being  the  boundary  layer  flow  over  a  swept  wing.  In  such  a  geometry  there 
is  a  flow  parallel  to  the  leading  edge  as  well  as  along  the  wings  chord.  In  flight  experi¬ 
ments  Gray  (1952)  noted  that  a  swept  wing  boundary  layer  reaches  a  transition  stage 
far  earlier  than  its  two-dimensional  counterpart.  The  explanation  of  this  was  given 
by  Gregory,  Stuart,  &:  Walker  (1955),  who  studied  the  archetypical  three-dimensional 
boundary  layer  above  a  rotating  disc.  By  considering  this  problem  from  both  a  the¬ 
oretical  and  experimental  standpoint  they  were  able  to  demonstrate  the  existence  of, 
what  are  now  known  as,  crossflow  vortices.  In  addition  to  crossflow  the  boundary  layer 
flow  over  a  swept  wing  also  brings  forward  the  question  of  how  an  applied  pressure 
gradient  will  affect  the  flows  stability. 

In  the  seminal  work  of  Taylor  (1923)  it  was  shown  that  centrifugal  forces  could 
support  a  toroidal  vortex  state  in  the  fluid  filled  gap  between  two  concentric  cylinders. 
Those  ideas  were  extended  by  Gortler  (1940)  and  applied  to  the  case  of  an  external 
boundary  layer  flow.  The  results  of  this  study  demonstrated  that  the  flow  over  a 
concave  plate  perpendicular  to  its  generator  is  susceptible  to  a  longitudinal  vortex 
state.  Such  instabilities  are  due  to  the  centrifugal  forces  present  in  situations  where 
wall  curvature  is  present.  Since  the  initial  work  of  Taylor  (1923)  and  Gortler  (1940) 
the  centrifugal  instability  problem  has  attracted  much  attention,  and  has  been  studied 
both  experimentally  and  theoretically,  using  mixture  of  asymptotic  and  numerical 
techniques. 

It  is  now  known  that  the  Gortler  vortex  instability  mechanism  is  a  truly  non-parallel 
phenomena,  in  that  it  is  not  possible  to  predict  a  growth  rate  locally  independent  of 
the  upstream  history  of  the  disturbance.  The  modes  initial  form  and  the  position  at 
which  it  is  imposed  both  play  a  significant  role  in  the  evolution  of  the  mode.  The  work 
of  Hall  (1983)  represents  the  first  attempt  at  solving  the  correct  equations  governing 
this  vortex  state.  A  Gortler  vortex  varies  over  commensurate  scales  in  the  wall  normal 
and  spanwise  directions,  but  evolves  on  the  same  scale  as  the  underlying  boundary 
layer  evolves  in  the  streamwise  coordinate.  Thus  in  the  limit  of  large  Reynolds  the 
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streamwise  derivative  of  the  pressure  and  the  streamwise  diffusion  terms  do  not  enter 
into  the  governing  equations.  However,  the  streamwise  derivatives  in  the  convective 
terms  cannot  be  neglected  as  they  are  the  same  order  of  magnitude  as  the  wall  normal 
and  spanwise  terms,  due  to  the  size  of  the  convective  velocities  in  these  directions. 
Thus  the  terms  that  can  legitimately  be  removed  render  the  system  parabolic  in  the 
streamwise  coordinate,  and  hence  enable  its  solution  by  a  relatively  cheap  marching 
procedure  from  some  given  streamwise  location  with  some  prescribed  initial  form.  This 
initial  form  must  be  consistent  with  the  governing  equations  which  impose  conditions 
on  the  normal  structure  of  the  disturbance;  these  conditions  are  derived  in  a  similar 
way  to  those  given  in  Goldstein’s  (1948)  study  of  a  flow  near  a  point  of  separation. 
The  parabolic  equations  can  then  be  marched  forward  and  the  magnitude  of  the  mode 
monitored;  when  the  streamwise  derivative  of  some  suitable  norm  of  the  mode  becomes 
zero  the  mode  is  said  to  be  at  its  neutral  location.  The  concave  curvature  of  the  wall 
is  allowed  to  vary  with  the  distance  from  the  leading  edge,  and  the  magnitude  of 
the  curvature  is  chosen  so  that  the  Gortler  mechanism  is  operable.  In  Hall  (1983) 
it  was  shown  that  as  Gortler  vortices  progress  downstream  in  order  for  them  to  be 
‘maintained’  the  wall  must  have  curvature  proportional  to  the  square  root  of  the 
distance  from  the  leading  edge.  It  was  observed  that  as  the  calculation  progresses 
the  growth  rate  curves  coalesce,  as  was  noted  and  exploited  in  Hall  (1982). 

In  an  asymptotic  study,  assuming  large  Gortler  number  (corresponding  to  the  be¬ 
haviour  far  downstream)  Hall  (1985)  showed  that  a  crossflow  will  stabilise  these  modes. 
In  that  study  far  larger  crossflows  can  be  described,  for  instance  Hall  (1985)  makes  pre- 
dictions  concerning  the  effect  of  increasing  the  crossflow  beyond  the  0(Re  '  )  bound, 
which  will  be  imposed  in  the  current  study.  (Here  Re  is  the,  large,  Reynolds  number  of 
the  flow).  In  Hall  (1985)  an  eigenvalue  problem  was  derived  which  allowed  the  Gortler 
number  required  for  neutrality  to  be  found.  As  the  degree  of  crossflow  increases  the 
Gortler  number  is  modified  and  the  vortex  activity  is  driven  to  a  thin  wall  layer,  until 
it  is  ultimately  extinguished.  This  behaviour  has  also  been  found  theoretically  in  the 
narrow  gap  Taylor  problem  by  Otto  &  Bassom  (1994a).  It  was  found  that  a  smaller 
crossflow  component  is  required  to  stabilise  the  most  unstable  Taylor  vortex  mode  than 
the  most  unstable  Gortler  modes,  so  it  was  conjectured  that  the  Gortler  mechanism 
would  be  more  prevalent  in  a  mixed  situation. 

In  a  preliminary  report  Otto  Denier  (1993a),  (henceforth  referred  to  as  OD), 
we  have  shown  that  crossflow  modifies  the  stability  characteristics  of  Gortler  vortices. 
In  that  study  several  assumptions  were  made  in  order  to  reduce  the  dimension  of 
the  parameter  space.  It  will  be  shown  that  the  results  obtained  in  OD  can  be  in 


2 


agreement  with  the  earlier  asymptotic  results  of  Hall  (1985),  depending  on  the  choice 
of  the  frequency  of  the  modes,  even  for  the  small  crossflows  considered  here. 

The  physical  problem  that  we  shall  consider  herein  is  the  flow  over  a  yawed  cylinder 
where  the  surface  of  the  body  possesses  a  degree  of  concave  curvature.  The  problem  is 
non- dimensionali sed  in  the  normal  manner,  using  the  dimensional  values  U0 0  (a  char¬ 
acteristic  free-stream  velocity),  L  (a  characteristic  lengthscale)  and  u  (the  kinematic 
viscosity).  There  are  two  key  non-dimensional  parameters, 

Re  =  G  =  2  R.U,  (1.1) 

the  former  being  the  Reynolds  number  which  is  assumed  to  be  large,  and  the  latter 
the  Gortler  number  which  is  assumed  to  be  an  order  one  quantity.  The  value  of  S,  a 
curvature  parameter,  is  chosen  to  tend  to  zero  as  Rt  — ►  oo  such  that  the  G  remains  of 
order  one,  this  ensures  that  the  Gortler  mechanism  is  operational.  For  a  comprehensive 
review  of  the  work  on  Gortler  vortices  the  reader  is  referred  to  Hall  (1990).  In  Bassom 
&  Hall  (1991)  a  study  was  made  of  the  connection  between  crossflow  and  Gortler 
vortices  within  the  high  Gortler  number  regimes.  The  problem  concerning  the  effect 
of  crossflow  on  the  Taylor  modes  was  discussed  by  DiPrima  &  Pridor  (1979). 

In  the  work  of  Hall  (1983)  the  growth  of  Gortler  vortices  within  a  Blasius  boundary 
layer  were  considered.  It  is  well  known  that  the  flow  over  a  flat  plate  is  self  similar  in 
the  downstream  coordinate,  x;  at  least  for  large  values  of  x  or  the  Reynolds  number. 
In  the  case  of  the  flow  over  a  swept  wing  there  is  a  natural  pressure  gradient  present 
and  it  is  for  this  reason  we  shall  consider  first  the  effect  of  this  pressure  gradient  on 
the  stability  characteristics  of  Gortler  vortices.  In  OD  it  was  shown  that  a  favourable 
pressure  gradient  promotes  the  instability  mechanism  and  is  was  conjectured  that 
this  is  connected  with  the  convective  nature  of  the  Gortler  modes.  We  will  assume 
that  the  pressure  gradient  takes  the  form  x2n,  so  that  in  the  free  stream  the  scaled 
streamwise  velocity  takes  the  form  xn.  This  form  of  the  pressure  gradient  leads  to 
the  Falkner-Skan  family  of  velocity  profiles  (Falkner  &  Skan  (1930)).  It  is,  of  course, 
possible  for  other  forms  of  pressure  gradient  to  be  used  but  these  would  only  serve 
to  complicate  the  problem  so  as  to  cloud  the  pertinent  issues.  In  the  work  of  Cooke 
(1950)  the  effect  of  yawing  on  this  problem  was  considered  and  it  was  found  that  the 
spanwise  momentum  equation  could  be  solved  after  obtaining  the  streamwise  profile. 
The  degree  of  crossflow  has  a  direct  correspondence  with  the  angle  of  sweep. 
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After  the  momentum  equations  have  been  linearized  about  the  basic  state  we  have 
four  coupled  partial  differential  equations  to  solve.  It  is  not  clear  how  to  impose  bound¬ 
ary  conditions  on  the  perturbation  pressure  in  this  primitive  formulation,  so  it  is  found 
to  be  convenient  to  eliminate  the  vortex  spanwise  velocity  and  pressure  components. 
We  first  Fourier  decompose  in  the  spanwise  direction  (where  the  spanwise  wavelength 
has  been  scaled  on  the  boundary  layer  thickness)  and  assume  that  the  modes  are  har¬ 
monic  in  the  temporal  variable  (we  note  that  it  is  beyond  the  scope  of  the  calculations 
to  be  presented  here  to  consider  modes  which  evolve  both  temporally  and  spatially). 
The  resulting  equations  then  form  a  coupled  system  of  partial  differential  equations  for 
the  velocity  components  as  functions  of  the  normal  and  streamwise  variables.  These 
equations  are  then  discretized  using  a  second  order  accurate  finite  difference  scheme 
in  the  normal  variable  and  a  Crank-Nicholson  scheme  in  the  streamwise  coordinate. 
A  stretched  grid  is  employed  in  the  normal  coordinate  so  as  to  retain  adequate  resolu¬ 
tion  near  the  wall  without  using  prohibitively  many  points  to  obtain  a  sufficiently  large 
outer  bound  for  the  calculations.  The  resulting  discretized  system  is  in  the  form  of  a 
coupled  penta-  and  tri-diagonal  system  which  was  inverted  directly  at  each  stream- 
wise  location  by  using  a  modified  form  of  a  Davis  Coupled  scheme.  This  scheme  has 
the  advantage  of  directly  inverting  the  system  (thus  yielding  a  true  implicit  scheme) 
without  incurring  the  extra  computational  costs  in  inverting  a  full  block  system.  This 
method  was  developed  for  a  weakly  nonlinear  calculation,  Bassom  &  Otto  (1993),  and 
has  been  shown  to  be  efficient  in  other  more  complex  calculations.  The  additional 
work  involved  in  implementing  the  Crank-Nicholson  scheme  for  the  streamwise  coor¬ 
dinate  was  minimal  and  did  lead  to  improved  accuracy  when  compared  to  the  results 
obtained  using  a  Backward  Euler  discretisation  in  the  streamwise  direction.  As  a  test, 
consistency  was  demonstrated  between  the  Euler  results  with  a  smaller  step-size  and 
the  Crank-Nicholson  results.  The  method  used  to  invert  the  system  is  totally  implicit; 
as  such  the  streamwise  step  is  restricted  to  be  Ay  rather  than  Ay 2  as  it  would  be 
in  an  explicit  scheme.  However  since  both  the  Euler  and  Crank-Nicholson  schemes 
axe  unconditionally  stable,  it  is  possible  to  use  a  larger  streamwise  steplength  without 
invalidating  our  results. 

Many  of  the  conclusions  found  here  have  also  been  observed  in  a  concurrent  study 
by  Zurigat  &  Malik  (1993).  That  work  concerns  steady  modes  in  situations  with 
favourable  pressure  gradients.  The  study  exploits  a  parallel  flow  approximation,  which 
although  erroneous,  shows  the  link  between  Gortler  modes  and  crossflow  modes  in 
three-dimensional  boundary  layers.  In  the  context  of  our  calculation  it  is  not  possible 
to  describe  the  crossflow  modes  (Gregory,  Stuaxt  &  Walker  (1955))  due  to  the  fact  that 
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the  streamwise  derivatives  that  have  been  removed  will  become  significant  in  such  a 
problem. 

The  paper  is  organized  as  follows:  in  section  2  we  derive  the  governing  equations 
and  discuss  the  conditions  used  to  derive  the  initial  conditions  for  the  calculation.  In 
section  3  we  discuss  the  numerical  techniques  used  to  advance  the  solution.  Section 
4  contains  a  summary  of  the  results,  section  5  includes  some  conclusions  and  possible 
future  extensions  to  this  work  and  in  an  appendix  we  present  some  asymptotic  results 
concerning  the  effect  of  pressure  gradient  on  the  Gortler  instability  mechanism. 

§2  Formulation 

We  consider  the  flow  over  a  yawed  cylinder  as  described  in  Hall  (1985)  in  which 
the  body  is  assumed  to  have  concave  curvature  perpendicular  to  its  generator.  Such 
a  flow  can  support  centrifugal  instabilities  and  we  wish  to  consider  the  effect  of  the 
angle  of  sweep,  which  has  a  direct  correspondence  to  the  amount  of  crossflow  present 
in  the  unperturbed  flow  field,  on  the  instability.  We  assume  that  the  pressure  gradient 
present  in  the  flow  is  of  a  form  which  is  conducive  to  a  self  similar  calculation,  namely 
x2n.  Since  p  ~  x2n,  the  scaled  velocity  at  the  extreme  edge  of  the  boundary  layer  is 
given  by  ue  ~  xn,  (from  the  balance  ututx  =  —  px).  The  viscous  boundary  layer  then 
grows  as  x1^ ,  and  we  introduce  the  similarity  variable  rj  =  yfx~*~ ,  where  y  is  the 
usual  scaled  boundary  layer  normal  variable.  Within  the  boundary  layer  the  flow  then 
takes  the  form 

( )  =  ^xn  f V,R7*  x^  » Re  79^j  ,  (2.1) 

where  the  functions  of  p,  f  and  g  satisfy, 

fnnn  +  =  n  (f*  -  l)  ,  gr,r,  +  ~fgn  =  0,  (2.2a,  b ) 

together  with  the  boundary  conditions 

/  =  /r?  =  9  =  0  at  ?7  =  0,  l,si  -►  A  as  p  ->  oo. 

If  A  =  0  these  equations  give  the  standard  Falkner-Skan  boundary  layer  profile.  In  the 
case  n  =  0  we  revert  to  the  Blasius  profile  in  which  case  g  is  a  linear  multiple  of  / 
and  the  Gortler  vortex  instability  mechanism  reverts  to  that  for  the  two  dimensional 
boundary  layer  flow  over  a  curved  plate  (the  reader  is  referred  to  Hall  (1985)  for  a 
discussion  of  this  case). 
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We  assume  that  the  bounding  surface  has  curvature  distribution  of  the  form  x(x)i 
where  \  is  a  positive  quantity  for  a  concave  wall.  Hall  (1982)  noted  that,  for  the  case 
of  a  Blasius  boundary  layer,  in  order  to  support  Gortler  vortices  far  downstream  of 
the  leading  edge  we  have  the  requirement  that  \  ~  x*-  I11  the  current  problem  this 
condition  is  now  modified  to  be  \  ~  x  ,  which  is  readily  deduced  from  the  balances 
required  near  the  right-hand  branch  of  the  neutral  curve. 

In  this  paper  we  restrict  our  attention  to  linear  disturbances  to  the  three-dimensioned 
boundary  layer.  We  perturb  the  basic  state  (2.1)  with  a  small  amplitude  wave  so  that 
the  total  flow  is  now  written  as 


q  =  (u,  Re  *tJ,  R7kw,p  +  R-'p)  +  A  (U,  R:>V,  Re  *  W,  Rj'P)  . 

Substituting  the  above  expansion  into  the  incompressible,  unsteady  Navier-Stokes 
equations  fitted  to  the  plate,  taking  first  the  limit  A  — ►  0  and  then  the  limit  Re  — ►  oo 
we  obtain  the  following  system  of  viscous  equations  governing  the  evolution  of  a  small 
amplitude  disturbance: 


~  ~  du  ~  du 

CU=UK+W 

(2.3a) 

=  U^  +  V^  +  GxaU+^, 
ax  ay  ay 

(2.36) 

-  ~  dw  ~  dw  dP 

cw  =  u—  +  v—  +  -f , 

dx  dy  dz 

(2.3c) 

dU  dV  dW 
dx  +  dy  +  dz 

(2.3  d) 

Here  the  differential  operator  C  is  given  by 

d2  &  _d  _d  _d  d 

C=  dy2+  dz 2  Udx~Vdy  Wdz~  dt' 

and,  anticipating  the  well  known  result  that  the  vortex  wavelength  scales  on  the  bound- 

1  /2 

ary  layer  thickness,  we  have  introduced  the  nondimensional  variable  z  =  Re'  z* / L. 
The  boundary  conditions  appropriate  to  (2.3)  are 


U  =  V  =  W  =  0  at  y  =  0,  oo, 


(2.4) 


corresponding  to  the  requirement  of  no-slip  at  the  wall  and  that  the  vortex  is  confined 
to  the  boundary  layer. 


At  this  juncture  we  note  there  is  no  streamwise  derivative  of  the  perturbation  pres¬ 
sure  and  no  streamwise  diffusion  in  the  perturbations  equations  (2.3)  thus  rendering 
the  system  parabolic  in  the  streamwise  coordinate  x.  As  such  the  system  (2.3)  can 
be  marched  downstream  with  some  suitable  initial  conditions  imposed  at  x  =  x  (say). 
However,  as  noted  by  Hall  (1983)  it  is  not  sufficient  to  choose  an  initial  form  for  the 
perturbation  velocity  fields  which  solely  satisfies  the  boundary  conditions  (2.4);  such 
a  choice  necessarily  introduces  singularities  at  x  =  x.  Indeed,  the  initial  perturbation 
quantities  must  satisfy  conditions  stronger  than  (2.4)  which  can  readily  be  derived 
by  following  the  method  outlined  in  Goldstein  (1948)  and  expanding  both  U  and  V 
in  Taylor  series  about  x  —  x.  Such  a  procedure  yields  the  result  that  the  initial 
perturbations  U  =  U(y)  and  V  =  V(y)  must  satisfy 

Uyyi 0)=0,  Uyyy  (0)  =  (a2  +  ift)  Uy  (0) ,  (0)  =  (2a2  +  »fl)  V„  (0) ,  (2.5) 

in  addition  to  the  boundary  conditions  (2.4).  Hall  (1983)  derives  the  condition  that 
U  must  satisfy  as  Y  — ►  oo,  which  when  modified  to  account  for  the  presence  of  the 
pressure  gradient  becomes 

6  ~ exp  r ')  ■  (26) 

As  noted  by  Hall  (1983)  a  parallel-flow  approximation,  in  which  the  term  udx  is 
neglected,  yields  the  large  y  structure 

U  ~  exp(-ay), 


thus  predicting  an  incorrect  decay  in  the  far  field  for  the  perturbation  quantities.  It 
is  for  precisely  this  reason  that  the  parallel  flow  approximation  is  erroneous  in  the 
context  of  the  Gortler  vortex  problem. 

In  order  to  reduce  the  system  (2.3)  to  a  form  which  is  more  convenient  from  a 
computational  standpoint  we  first  Fourier  decompose  in  the  spanwise  and  temporal 
variable  such  that  dz  — >  ia  and  <9<  — ►  —ifl,  where  both  a  and  O  are  real  quantities, 
and  eliminate  the  perturbation  pressure  and  spanwise  velocity  perturbation  to  give 


(£-) 


d2U  du  dU  d2u 


V  +  Ga  xuU  +  2  — — b  2 


dxdy  dx  dx  dxdy 


-  5s fZ  dVd2u  dV  d2u 

.1.  ij - 1 - 1 - 

dx2dy  dx  dy2  dy  dxdy 

-  d?u  dvd2U  .  d2w  ~  dwdU 
+  V dxdy2  +  dx  dy2  +2tadxdyU+  ^ dx  dy 

dv  d2V  .  ~  d2  w  2  tt  dv  2  ~  dv 

-dfW+,aVW  +  aUd;  +  aVd;=0’ 
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together  with  equation  (2.3a),  where  the  differential  operator  £  is  now  given  by 


£  = 


£_ 

dy 2 


a  —  u— - v— - taw  +  ill. 


dx 


dy 


The  problem  now  reduces  to  the  numerical  solution  of  the  governing  equations  (2.3a) 
and  (2.7)  subject  to  the  boundary  conditions 


dV 

U  =  V  =  — —  =  0  at  y  =  0, 
dy 

~  ~  dV 

U,Vy— - >0  as  y  — >  oo. 

dy 


(2.7) 


given  some  initial  perturbation  at  i  =  i,  consistent  with  the  conditions  (2.5)  and 
(2.6),  marching  downstream  in  x.  In  the  next  section  we  present  such  a  numerical 
procedure. 


§3  Numerical  Methods 

The  method  used  to  numerically  integrate  the  governing  equations  (2.3a)  and  (2.7) 
parallels  that  used  by  Hall  (1983),  however  as  there  are  some  significant  differences  we 
now  give  an  indication  of  the  pertinent  points  regarding  our  numerical  scheme. 

Following  Hall  (1983)  we  elect  to  solve  the  system  in  terms  of  the  similarity  variable 
T)  =  y/x  2  .  We  then  write  U  =  U(x,  rj),  together  with  similar  expressions  for  the  ver¬ 
tical  component  of  the  perturbation  velocity  field,  and  take  our  underlying  boundary 
layer  flow  to  be  described  by  (2.2a, b).  The  system  can  then  be  transformed  from  the 
(x,y)-space  to  the  (x,y)- space  by  using  the  fact  that 

d  d  1  +  nrj  d  d  Id 

dx  dx  2  x  dr/'  dy  x dr] ’ 

the  higher  order  derivatives  transform  in  a  similar  way  as 

q2  _  j;  o2  q3  __  _  <a3 

uyy  ~  -i -nuvv  °yyy  ~  3-3n  uiynni 

%  X  2 

dly  =  dlx  =  ~  U  ~hl  (d1  +  vtfv)  +  -Eltfz, 

&  X  2  X  2 

%yy  =  ( 2d In  + 

Such  a  transformation  allows  the  calculation  to  spread  in  the  normal  extent  with  the 
growth  of  the  boundary  layer.  In  order  to  achieve  adequate  resolution  at  the  wall, 
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whilst  still  retaining  a  sufficiently  large  outer  bound  for  the  calculations  and  without 
using  a  prohibitively  large  number  of  normal  grid-points,  an  algebraic  stretching  of 
the  form  introduced  by  Macaraeg,  Streett  &  Hussaini  (1988)  was  used  for  the  r?  grid. 
In  the  majority  of  the  calculations  presented  herein  150  points  were  used  in  the  normal 
coordinate,  with  an  outer  bound  at  77  =  20.  The  basic  flow  quantities  u,  v  and  w  were 
obtained  a  priori  by  using  a  standard  fourth  order  Runge-Kutta  technique,  coupled  to 
a  two-dimensional  real  secant  method,  to  integrate  the  similar  equations  (2.2);  these 
were  then  projected  onto  the  similarity  grid  at  each  downstream  station. 

In  order  to  march  the  equations  downstream  from  an  initial  position  x  =  x  a, 
formally  second  order  accurate,  Crank- Nicholson  scheme  was  used  in  the  downstream 
coordinate.  If  we  compare  this  scheme  with  a  Back  war  d-Euler  scheme,  this  being  the 
method  used  by  Hall  (1983)  and  others,  we  find  that  in  using  the  two  schemes  for 
ff  +  =  0,  we  are  required  to  solve 

Backward-  Euler: 

<f>+  dP<f>+  _  <t>~ 
e  +  dyP  t  ’ 

Crank-Nicholson: 

1  dP<f>+  <f>~  1  dP<f>~ 

e  +  2  dyP  ~  e  2  dyP  ' 

for  <£+  as  a  function  of  y.  Here  the  superscripts  +  and  —  denote  evaluation  at  the 
current  and  previous  x-position  respectively,  and  e  is  the  streamwise  step-length.  We 
note  that  although  the  system  which  requires  inversion  to  determine  <f>+  is  of  the  same 
form  in  both  case,  the  extra  accuracy  comes  from  evaluating  the  at  the  mid-point 
and  hence  the  streamwise  derivative  then  becomes  a  centred  difference.  The  only 
additional  computation  that  results  is  involved  in  computing  the  term  %£,•  which  is 
outweighed  by  the  additional  accuracy  obtained.  In  all  the  calculations  reported  here 
a  streamwise  step-length  of  0.1  was  used. 

The  system  is  now  discretized  using  a  second  order  accurate  (for  the  fourth  deriva¬ 
tive,  this  requires  five  points)  finite  difference  scheme  in  the  normal  7?  coordinate  which 
when  combined  with  the  Crank-Nicholson  discretisation  for  the  streamwise  derivatives 
allows  the  system  to  be  written  in  the  form 

anVn+2+bnVn+i+CnVn+  dnVn^+enVn-2+fnUn+l+gnUn  +  hnUn-i  =  R^\  (3.1a) 
and 

®n^n+l  d"  bnUn  -j-  CnUn_i  -|-  <2-n^n  =  ^ •  (3.16) 
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Here  the  index  n  denotes  evaluation  at  the  nth  normal  grid-point,  the  terms  on  the 
left-hand-sides  of  these  expression  are  evaluated  at  x  =  x  +  e  and  the  and 
denotes  terms  evaluated  at  the  previous  streamwise  position.  (For  the  sake  of  brevity 
the  form  of  the  coefficients  in  (3.1a, b)  Eire  not  presented  here;  they  can  be  obtained  on 
application  to  the  authors). 

Here  we  highlight  the  major  difference  between  the  present  scheme  and  that  used  in 
Hall  (1983)  (apart  from  the  additional  accuracy  obtained  in  the  present  calculations 
due  to  the  Crank-Nicholson  discretisation  of  the  streamwise  derivatives).  The  method 
used  in  Hall  (1983)  is  to  move  the  terms  in  (3.1a)  associated  with  u  to  the  right  hand 
side,  hence  evaluating  them  at  the  previous  point,  and  then  invert  the  penta-diagonal 
system  for  v.  Using  this  solution  it  is  then  possible  to  solve  the  tri-diagonal  system 
(3.1b)  for  u.  Here  we  choose  to  invert  the  entire  system  in  one  calculation,  thus 
preserving  the  structure  of  the  system  and  rendering  our  scheme  truly  second  order 
accurate  in  the  streamwise  coordinate.  The  method  used  for  this  inversion  represents 
a  modified  form  of  a  Davis  Coupled  Scheme  and  has  been  outlined  in  Otto  &  Bassom 
(1993)  (the  reader  is  referred  to  that  paper  for  full  details). 

To  monitor  the  growth  (or  decay)  of  an  initial  disturbance  as  the  calculation  pro¬ 
ceeds  downstream  we  define  the  energy  as 


If  we  define  the  growth  rate  as 


Eiv)  1  -  n 
EM  +  2x  ’ 


(3.2) 


then  a  neutral  point  is  the  x-station  where  the  real  part  of  the  quantity  cr,  defined  in 
(3.2),  changes  sign.  The  imaginary  part  of  a  gives  us  the  wavespeed  in  the  streamwise 
coordinate.  It  should  be  noted  that  for  the  steady  Gortler  problem,  ft  =  0,  we  have 
Imag  (cr)  =  0. 


§4  Results  and  Discussion 

In  a  preliminary  report  (Otto  h  Denier  (1993a))  we  presented  some  limited  re¬ 
sults  concerning  the  effect  of  crossflow  on  the  Gortler  instability  mechanism  at  0(1) 
wavelengths;  however  several  choices  of  various  parameters  were  made  which  limits  the 
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applicability  of  that  study.  Here  we  shall  now  endeavour  to  provide  a  fuller  description 
of  the  effect  of  crossflow  on  the  Gortler  vortex  instability  mechanism. 

In  order  to  decrease  the  number  of  parameters  in  our  study  we  will,  throughout  the 
calculations  reported,  fix  the  Gortler  number  and  note  that  for  larger  Gortler  numbers 
the  results  will  be  quantitatively  similar  to  those  presented  here  (keeping  in  mind 
that  for  larger  Gortler  numbers  the  flow  becomes  unstable  at  a  smaller  value  of  the 
streamwise  variable  x).  The  value  chosen  for  the  Gortler  number  is  G  =  0.05,  this 
value  is  taken  solely  for  convenience.  The  work  of  Hall  (1983)  clearly  demonstrates, 
in  the  case  of  Gortler  vortices  in  two-dimensional  boundary  layers,  that  the  initial 
form  of  the  perturbation,  and  the  position  at  which  it  is  imposed,  are  key  factors  in 
determining  the  stability  characteristics  of  flow  field.  Not  surprisingly  such  is  the  case 
in  the  present  problem,  however  as  we  are  predominantly  concerned  with  the  effect  of 
both  an  imposed  pressure  gradient  and  crossflow  we  choose  to  ignore  this  effect  and 
consider  only  the  case  when  the  initial  perturbation  is  of  the  form 

(U,V)  =  o), 

(note  that  this  is  precisely  the  form  used  by  Hall  (1983)).  Here  T?ext  is  a  measure  of 
the  normal  extent  of  the  initial  perturbation;  for  the  sake  of  definiteness  we  choose 
»?ext  =  4  for  all  the  calculations  reported  here.  This  condition  is  now  imposed  at  x  =  20 
and  the  numerical  scheme  described  above  is  used  to  march  the  solution  in  x  at  each 
stage  monitoring  the  growth  rate  a  of  the  disturbance.  We  note  in  passing  that  the 
value  at  which  the  initial  perturbation  is  imposed  was  varied  ( x  =  30  and  x  =  40); 
in  all  case  the  growth  rate  curves  coalesce  into  the  unique  right-hand  branch  of  the 
neutral  curve  far  downstream.  It  is  in  this  regime  that  the  parallel  flow  approximation 
achieves  validity,  however  as  demonstrated  by  Hall  (1982)  this  asymptotic  regime  can 
readily  be  described  without  recourse  to  this  ad  hoc  approximation  technique. 

One  of  the  major  differences  between  this  work  and  OD  is  that  here  we  shall  retain 
the  variation  of  x  with  x.  Thus  for  each  pressure  gradient  two  forms  of  wall  curvature 
will  be  considered;  (:)  x  =  1  (OD),  and  (it)  x  —  x~*~ /x~^  (this  being  the  wall 
curvature  required  to  support  Gortler  vortices  far  downstream). 

We  consider  first  the  effect  that  the  modification  of  the  basic  flow  due  to  the  presence 
of  the  pressure  gradient  has  upon  the  stability  characteristics  in  the  absence  of  the 
crossflow  component  of  the  basic  flow.  As  was  demonstrated  by  OD,  in  this  case  the 
most  unstable  mode  is  found  to  be  steady;  we  therefore  set  =  0.  In  figure  1  we  show 
three  neutral  curves  for  the  cases  n  =  1/20,  n  =  0  and  n  =  —1/20  for  the  choice  (i)  of 
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constant  curvature;  these  correspond  to  a  favourable-,  no-,  and  an  adverse-pressure 
gradients  respectively.  The  results  are  presented  in  terms  of  the  local  Gortler  number 
and  local  wavenumber  which  are  defined  as 

Gx  —  Gx^,  ax  =  ax  *  ,  (4.1) 

and  have  been  derived  with  respect  to  the  local  boundary  layer  thickness.  The  neutral 
curves  were  generated  by  fixing  the  Gortler  number  G  =  0.05  (see  the  discussion 
above),  choosing  a  spanwise  wavenumber  a  and  marching  the  solution  to  (2. 3a, 2. 7) 
downstream  until  such  point  at  which  the  growth  rate  a  (given  by  (3.2))  changes  sign. 
At  this  point  the  local  Gortler  number  and  local  wavenumber  can  be  calculated  from 
(4.1).  This  procedure  is  then  repeated  for  different  values  of  spanwise  wavenumber  a 
to  produce  a  curve  of  the  form  given  in  Figure  1. 

In  figure  2  we  present  the  neutral  curves  for  the  curvature  distribution  (ii),  again 
with  G  =  0.05.  In  both  cases  it  is  readily  seen  that  for  the  case  of  a  favourable  pressure 
gradient  the  range  of  local  wavenumber  over  which  the  flow  is  unstable  is  increased 
and,  additionally,  the  minimum  value  of  the  local  Gortler  number  Gx  is  decreased 
with  increasing  n.  We  note  that  this  thickening  of  the  neutral  curve  corresponds  to 
an  increase  in  the  streamwise  extent  over  which  a  disturbance  of  given  wavenumber 
a  is  unstable.  Thus,  for  the  two  representative  curvature  distributions  considered 
here,  it  is  the  case  of  boundary  layer  with  a  favourable  pressure  gradient  that  is  most 
susceptible  to  Gortler  vortices.  In  the  limit  of  large  spanwise  wavenumber  a  and  large 
Gortler  number  G  it  is  possible  to  describe,  in  a  self  consistent  asymptotic  setting,  the 
structure  of  the  Gortler  vortex  in  the  vicinity  of  the  right-hand  branch  of  the  neutral 
curve.  The  details  of  the  calculation  relevant  to  the  right-hand-branch  of  the  neutral 
curve  are  given  in  the  appendix  at  the  end  of  this  paper. 

In  figure  3  we  present  plots  of  the  (scaled)  growth  rate,  fl* ,  and  the  position  of 
vortex  activity,  77*,  as  a  function  of  the  scaled  wavenumber  <7*  for  various  values  of  n 
(see  the  appendix  for  the  precise  definition  of  these  scaled  quantities).  From  figure  3 
we  note,  as  found  by  Denier  et  al  (1991),  that  /?*  ~  l/<r*2  as  a*  — ►  0;  this  fact  can 
then  be  used  to  determine  the  most  unstable  Gortler  vortex  mode  as  in  Denier  et  al. 
To  clearly  see  the  effect  of  changing  the  pressure  gradient  (by  varying  the  parameter 
n )  on  the  right-hand-branch  of  the  neutral  curve  we  present,  in  figure  4,  a  plot  of 
the  scaled  neutral  wavenumber  a*,  at  which  ft*  =  0,  together  with  the  position  of 
vortex  activity  77*.  As  n  decreases  the  scaled  wavenumber  <7*  decreases;  thus  a  less 
favourable  pressure  gradient  moves  the  right-hand-branch  of  the  neutral  curve  to  the 
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left.  In  this  sense  a  boundary  layer  flow  with  a  favourable  pressure  is  more  unstable 
to  Gortler  vortices.  Similarly,  for  decreasing  n,  the  position  r?*  at  which  the  neutral 
vortex  resides  increases.  Combining  the  results  of  our  numerical  calculations  with  the 
large  wavenumber  asymptotics  which  appear  in  the  appendix  we  can  clearly  see  that 
the  range  of  local  wavenumbers  over  which  the  disturbance  is  growing  is  increased  and 
also  the  critical  local  Gortler  number  is  decreased.  This  result  is  in  stark  contrast 
to  the  case  of  Tollmien-Schlichting  waves  in  an  external  boundary  layer  where  it  has 
been  demonstrated  that  an  adverse  pressure  gradient  is  destabilizing  (see  Rosenhead 
(1963),  Ch.  IX).  We  conjecture  that  such  a  destabilization  due  to  a  favourable  pressure 
gradient  is  due  to  the  convective  nature  of  the  Gortler  vortex  instability.  However, 
since  the  question  of  whether  the  Gortler  vortex  instability  is  convective  or  absolute 
is  still  undecided,  although  experimental  observations  do  support  the  commonly  held 
belief  that  it  is  in  fact  a  convective  instability,  we  believe  that  this  point  warrants 
further  investigation,  the  results  of  which  we  hope  to  report  in  a  forthcoming  article. 

We  note,  from  figure  1  and  2,  the  apparent  absence  of  a  left-hand  branch  of  the 
neutral  curves.  This  arises  due  to  the  computational  difficulty  in  obtain  neutral  values 
relevant  to  this  branch  due  to  the  fact  that  exceptionally  small  es  of  the  wavenum¬ 
ber  a  are  required  in  such  a  calculation.  At  this  point  we  emphasize  that  the  evolution 
of  a  single  mode  of  fixed  spanwise  wavenumber  a  the  local  wavenumber  varies  accord¬ 
ing  to  (4.1)  as  the  calculation  proceeds  downstream.  In  this  case  the  loci  of  local 
wavenumber  is  given  by  the  dashed  line  in  figure  2  and  as  such  there  is  a  finite  interval 
in  the  streamwise  direction  for  which  oy  is  positive  and  the  mode  is  growing  in  am¬ 
plitude.  In  this  interval  we  would  then  expect  nonlinear  terms  to  come  into  play  and 
the  results  of  Hall  (1988),  Otto  &  Bassom  (1994b)  and  Denier  &  Hall  (1993)  become 
relevant.  The  loci  shown  in  figure  2  is  for  the  representative  case  a  =  0.185;  the  solid 
block  shows  the  position  at  which  cry  first  changes  sign  and  we  note  that  this  trajectory 
does  not  cross  the  right  hand  branch  with  the  range  of  the  current  calculation  (4000 
steps  of  length  0.1  in  x). 

We  now  turn  our  attention  to  the  effect  of  cross  flow  on  the  Gortler  instability.  In 
the  asymptotic  regime  relevant  to  the  right-hand  branch  of  the  neutral  curve  (ie  large 
wavenumber  a  and  large  Gortler  number  G )  Hall  (1985)  demonstrated  that  the  effect 
of  crossflow  is  to  shift  the  right-hand  branch  of  the  neutral  curve  to  the  left.  We  will 
show  that  such  an  effect  can  be  achieved  even  with  the  modest  values  of  crossflow. 

In  our  earlier  study,  Otto  &  Denier  (1993a),  the  effect  of  crossflow  was  considered 
for  the  three  forms  of  pressure  gradients  given  above  and  it  was  shown  that  all  cases 
were  affected  in  a  similar  manner.  For  this  reason  we  choose  to  consider  only  the 
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case  of  the  favourable  pressure  gradient  here  as  this  is  the  most  unstable  in  the  sense 
described  above.  In  figure  5  we  show  a  different  form  of  neutral  curve,  where  we  fix 
the  wavenumber  at  a  =  0.195  (with  the  other  parameters  as  n  =  1/20,  G  =  1/5  and 
X  =  1)-  These  choices  were  made  so  that  the  calculation  is  in  a  parameter  regime 
which  is  relevant  to  the  right-hand  branch  of  the  neutral  curve.  In  this  case  the 
“stability  bubbles” are  actually  closed.  Here  we  present  results  for  increasing  values 
of  the  crossflow  parameter  A,  the  curves  shown  are  for  A  =  0,0.1, 0.2.  The  triangular 
symbols  depict  the  points  at  which  the  disturbance  begins  to  grow,  oT  first  becomes 
positive,  and  the  circular  symbols  where  it  first  begins  to  decay,  ar  becomes  negative. 
We  notice  that  the  effect  of  increasing  the  crossflow  parameter  A  is  to  shift  these  '  ves 
rather  than  compressing  them  as  might  be  expected  from  the  earlier  asymptotic  >s 

of  Hall  (1985)  and  Otto  &:  Bassom  (1993). 

In  order  to  explain  this  apparent  discrepancy  we  revert  to  the  other  form  of  neutral 
curve,  that  is  where  we  fix  Q  and  vary  the  wavenumber.  One  of  the  difficulties  in 
comparing  our  results  with  the  asymptotic  predictions  of  previous  authors  is  that  we 
must,  necessarily,  fix  the  frequency  of  the  disturbance  at  the  outset  for  our  numerical 
calculation.  In  figure  6  we  show  the  calculation  for  Q  —  1/20  with  A  =  1  and  Q  =  A  =  0 
for  G  =  1/20,  x  =  1  and  n  =  1/20.  In  this  case  the  right-hand  branch  for  the  case  with 
crossflow  has  been  moved  to  the  left,  and  now  lies  inside  the  neutral  curve  for  the  two- 
dimensional  boundary  layer  case.  However,  we  note  that  the  critical  neutral  Gortler 
number  is  only  slightly  modified  by  the  presence  of  the  crossflow.  Thus,  although 
crossflow  has  an  appreciable  effect  of  the  right-hand  branch  structure,  a  result  that  is 
apparent  from  the  asymptotic  results  of  Hall  (1985)  and  Otto  &  Bassom  (1993),  it  does 
not  appear,  on  the  basis  of  our  calculations,  to  completely  destroy  unstable  Gortler 
vortices.  In  particult -,  for  a  given  flow  situation  conclusions  concerning  the  instability 
of  a  curved  boundary  layer  in  the  presence  of  crossflow  can  only  be  determined  by 
solving  the  full  governing  equations  (2.3)  starting  from  some  suitable  initial  condition. 

To  complete  our  discussion,  in  figure  7,  we  present  the  streamwise  variation  of  the 
growth  rate  cr,  for  the  particular  choice  of  parameters  A  =  0.2,  a  =  0.195,  G  =  0.05, 
n  =  1/20  and  x  =  1>  from  which  we  can  readily  see  that  the  growth  rate  ar  is  positive 
over  a  large  streamwise  interval.  Finally  in  figure  8  we  show  the  contours  of  the  real 
part  of  the  normal  perturbation  velocity  component  V  at  z  =  0,  and  in  figure  9  we 
shown  the  contours  of  the  real  part  of  V  at  r)  =  1.  Figure  8  clearly  demonstrates  the 
growth  of  the  vortex  motion  as  the  calculation  proceeds  downstream  and  also  highlights 
the  presence  of  the  imaginary  part  of  the  growth  rate  a  on  the  vortex  motion  while 
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from  figure  9  we  see  the  vortices  migrate  in  the  x-z  plane  as  expected  in  the  presence 
of  crossflow. 

§5  Conclusions 

In  this  study  we  considered  the  effect  of  crossflow  on  the  Gortler  vortex  instability 
mechanism.  Perhaps  just  as  important  as  the  conclusions  concerning  crossflow,  is  the 
finding  that  the  favourable  pressure  gradient  flows  lead  to  far  more  unstable  Gortler 
modes,  a  fact  that  we  contribute  to  the  convective  nature  of  the  Gortler  instabilities. 
We  have  limited  ourselves  to  calculations  in  which  the  crossflow  is  small,  this  is  nec¬ 
essary  so  that  waves  with  spanwise  wavelength  comparable  with  the  boundary  layer 
thickness  can  be  considered.  It  is  found  that  the  dominant  effect  of  crossflow  is  to 
re-orientate  the  modes  and  change  the  streamwise  extent  over  which  they  grow.  In 
the  conventional  Gortler  problem  it  is  well  known  that  the  modes  are  steady,  however 
in  this  three-dimensional  case  the  ‘most  unstable’  mode  is  not  necessarily  steady.  We 
note  that  as  the  crossflow  is  further  increased  the  modes  start  to  grow  after  shorter 
streamwise  distances,  until  the  terms  d\  and  dP/dx  become  significant.  The  modes 
characteristics  can  then  be  determined  only  by  solving  the  full  linearized  equations.  In 
this  case  these  modes  can  be  thought  of  as  crossflow  modes  modified  by  curvature.  It  is 
now  accepted  that  good  approximations  to  these  equations  can  be  found  using  a  form 
of  the,  so  called,  Parabolised  Stability  Equations  (PSE),  see  for  example  Bertolotti 
(1991).  It  is  worth  emphasising  that  the  PSE  are  exact  for  the  two-dimensional  Gortler 
problem,  since  in  this  case  the  streamwise  wavenumber  is  zero. 

It  would  be  interesting  to  determine  how  the  crossflow  will  affect  the  receptivity 
of  a  situation  to  Gortler  vortices.  We  presume  that  the  ‘extra’  shear  caused  by  the 
presence  of  the  crossflow  will  increase  the  receptivity  coefficient.  In  a  recent  review 
article  by  Bassom  &  Seddougui  (1994),  the  receptivity  coefficient,  which  is  a  measure  of 
the  coupling  between  an  external  perturbation  (for  example,  in  the  form  of  an  isolated 
surface  roughness  element)  and  the  resulting  flow  disturbance,  was  calculated  for  the 
most  unstable  Gortler  vortices  and  it  was  found  to  increase  with  increasing  crossflow. 
Such  a  receptivity  calculation  would  in  a  similar  fashion  to  that  presented  by  Denier, 
Hall  &  Seddougui  (1991),  using  the  modifications  to  the  initial  conditions  recently 
proposed  by  Bassom  &  Hall  (1994). 

Another  interesting  question  concerns  the  effect  that  the  crossflow  will  have  on  the 
secondary  stability  of  a  boundary  layer  which  has  been  perturbed  by  the  presence  of  a 
Gortler  vortex.  In  Hall  &  Horseman  (1991)  and  Otto  &:  Denier  (1993b),  the  temporal 
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inviscid  stability  of  two  nonlinear  vortex  states  were  considered,  derived  in  Hall  (1988) 
and  Otto  &  Bassom  (1994b)  respectively.  In  the  present  calculation  we  have  only  a 
linear  vortex  state  however  the  stability  of  the  boundary  layer  with  this  imposed  vortex 
motion  can  still  be  analysed.  In  the  case  of  a  Gortler  vortex  in  the  presence  of  a  small 
crossflow  the  problem  of  the  secondary  instability  can  be  posed  in  am  asymptoticadly 
rigorous  setting  amd  the  resulting  three-dimensional  Rayleigh  equation  can  be  solved 
using  the  methods  developed  in  Otto  &  Denier  (1993b).  The  results  of  the  study 
of  Otto  &  Denier  (1993b)  have  demonstrated  the  existence  of  three  inviscid  modes 
which  cam  lead  to  the  flow’s  breakdown.  Firstly  there  are  the  even  amd  odd  modes, 
which  give  rise  to  the  horseshoe  and  wavy  vortices  respectively.  There  is  adso  a  second 
even  mode  which  has  a  larger  temporad  growth  rate  amd  persists  over  a  larger  ramge  of 
secondary  streaimwise  wavelengths.  This  mode  has  a  large  second  haurmonic  component 
the  magnitude  of  which  is  actually  equal  to  that  of  the  first  haurmonic  at  the  waill.  A 
similar  mode  has  adso  been  observed  in  stability  of  the  flow  in  a  comer  flow  by  Dhamadc 
(1993).  Studies  of  both  the  receptivity  and  secondary  instability  problem  are  currently 
underway  and  we  hope  to  report  on  the  results  of  these  studies  in  the  neaur  future. 
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Appendix 

Here  we  present  some  details  concerning  the  effect  of  the  pressure  gradient  on  the 
Gortler  instability  mechanism.  Noting  that  the  results  presented  in  §3  axe  relevant  to 
the  case  when  both  the  wavenumber  a  and  the  Gortler  number  G  are  0(1)  quantities 
we  will  restrict  our  attention  to  the  case  of  the  large  Gortler  number  limit. 

To  proceed  we  consider  first  the  asymptotic  regime  appropriate  to  the  right-hand- 
branch  of  the  neutral  curve  (we  will  follow  closely  the  work  of  Denier  et  al  (1991);  the 
reader  is  referred  to  that  paper  for  full  details).  Thus,  setting  A  =  $1  —  0  in  (2. 3a, 2. 7) 
we  seek  solutions  in  the  limit  G  — ►  oo  with  a  =  0(G1^4).  Thus  writing 

a  =  <rG1/4  (A.l) 

we  seek  a  solution  of  (2. 3a, 7)  confined  to  a  layer  of  depth  fc-1/2  ~  G'-1/8  centred  on 
the  location  y  =  y(x).  The  appropriate  expansion  is  given  by 

tt  =  |u0(x,O  +  fc~1/2“i(aJ»O  +  '”}exP  [Gl/2y  P(x)dx  » 

where  £  =  <r1/2G1/8(y  —  y).  The  eigenrelation  is  then  found  to  be  given  by 

[(«P  +  a2)2  ~  x«  %]  |,=9  =  0,  (A. 2) 

together  with  the  condition  that  ensures  the  vortex  is  confined  to  the  layer  centred  on 

y  =  y; 

[(up  +  (T2)2  -  X««»]  I y=9  =  0.  (A. 3) 

Thus  the  conditions  (A. 2)  and  (A. 3)  serve  to  determine  both  the  growth  rate  /3  and 
the  position  y  as  a  function  of  the  seeded  wavenumber  a.  For  the  self-similar  boundary 
layer  (2.2a)  we  have  u  =  xn/'(r/),  where  rj  =  y/x  2".  Writing 

P  =  /?*  (x*V )  1  ,  °  =  *n/ 2  (A.4) 

the  eigenrelation  (A. 2)  together  with  (A. 3)  can  be  rewritten  as 

(/'r + v*? = wm'D' + <e2)  =  (/’/")'  (a.5) 

where  /'  is  to  be  evaluated  at  some  position  77*  where  (A. 5b)  holds.  We  note  that 
n  does  not  appear  explicitly  in  (A. 5);  its  presence  is  felt  through  the  effect  on  the 
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boundary  later  profile  /  (see  (2.2a)).  Plots  of  /3*  versus  a*  are  given  in  the  text  of  the 
paper. 

Here  we  note  that  there  a  a  critical  value  of  n  =  nc  <  0  at  which  the  boundary 
layer  equations  (2.2)  no  longer  yield  physically  relevant  solutions;  indeed  at  this  point 
we  have  /"( 0)  — ►  0  as  n  — ►  nc.  For  n  <  nc  the  boundary  layer  equations  then  predict 
solutions  which  have  a  finite  region  of  separated  flow  for  all  values  of  the  streamwise 
coordinate  a:;  such  a  solution  is  of  no  physical  relevance. 

To  complete  the  description  of  the  wavenumber  spectrum  for  Gortler  vortices  in  the 
large  Gortler  number  limit  we  can  follow  Denier  et  al  (1991)  and  consider  the  inviscid 
Gortler  vortex  modes  which  lie  in  the  parameter  regime  G  1  with  a  =  0(1).  In  this 
case  there  is  an  exact  solution  of  the  linearized  governing  equations  (see  Denier  et  al 
(1991)  for  full  details)  in  which  the  growth  rate  G_1/2/3  is  given  by 

G-’/^  =  f, 

which  depends  only  on  the  curvature  and  the  spanwise  wavenumber.  Thus  the  growth 
rate  of  the  inviscid  Gortler  vortex  is  unaffected  by  the  presence  of  a  pressure  gradient; 
note  however,  that  the  normal  structure  of  the  inviscid  mode  is  affected  due  to  its 
dependence  on  the  basic  flow  u  (see  Denier  et  al  (1991)).  It  is  then  the  inviscid  regime 
which  represents  a  “turning  point”  with  respect  to  the  effect  of  pressure  gradient 
on  the  Gortler  vortex  instability.  For  a  1  (with  G  1)  a  favourable  pressure 
gradient  destabilizes  the  flow  whereas  for  a  <C  1  (again  with  G  >■  1)  a  favourable 
pressure  gradient  stabilizes  the  flow.  This  last  point  can  be  derived  from  the  small 
wavenumber  asympotitics  presented  in  Hall  &  Morris  (1992)  and  is  left  as  an  exercise 
for  the  interested  reader. 
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Figure  1:  Neutral  curves  for  case  (i),  with  n  =  1/20,  0  and  n  =  —1/20,  they  curves  appear 
in  ascending  order,  that  is  the  highest  curve  corresponds  to  n  =  —1/20. 


Figure  2:  Neutral  curves  for  case  (*»),  with  n  =  1/20,  0  and  n  =  —1/20,  they  curves 
appear  in  ascending  order. 
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Figure  8:  Contours  of  the  real  part  of  V  at  z  =  0. 
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